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Abstract 

A point particle treatment to the statistical mechanics of BPS black holes in 
Einstein-Maxwell-dilaton theory is developed. Because of the absence of the 
static potential, the canonical partition function for N BPS black holes can 
be expressed by the volume of the moduli space for them. We estimate the 
equation of state for a classical gas of BPS black holes by Pade approxima- 
tion and find that the result agrees with the one obtained by the mean-field 
approximation . 
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I. INTRODUCTION 



The thermodynamics of self-gravitating system has some importance in the context of 
cosmology. The long-range nature of gravity gives rise to the breakdown of the conven- 
tional description of statistical properties, such as the additivity of the free energy or other 
'extensive' thermodynamical quantities. While the statistical treatment of non-equilibrium 
systems has recently been developed, the existence of the long range force seems essential 
for explaining the formation of the peculiar local structure in the system. 

Recently, de Vega and Sanchez studied the statistical mechanics of the self-gravitating 
particle gas Q. They showed that the 'thermodynamic limit' must be taken as N — > oo with 
JV/V 1 ' 3 fixed in three dimensions, where N is the number of particles and V is the volume 
of the system. This treatment is required by the long-range nature of gravity, because there 
is no ad hoc cut-off scale. [] 

On the other hand, two of the present authors have studied the statistical mechanics of 
the well-separated charged particles in Einstein-Maxwell-scalar theory [fj]. The system was 
found to be unstable in general, though there exists a critical case in which the static forces 
are cancelled each other. 

In general relativity, the critical case was investigated and exact static solutions which 
describe multi-black hole systems have been obtained |§-|S|]. In such a case, the attractive 
forces (the gravitational and scalar-mediated forces) and the repulsive forces (the electro- 
static force) between black holes are exactly cancelled in the static limit. In this case, the 
energy of the system was calculated for small velocities but for small distances between the 
black holes 0-0. 

In this paper, we investigate the multi-black hole system by adopting the technique of 
de Vega and Sanchez IJ. Note that a mass-charge relation is satisfied for each individual 



1 Actually, the horizon length is a candidate for the cut-off scale in the universe. 
2 For various specific models and their applications, see [13 23]. 
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black hole. Such extreme black holes, or so-called BPS black holes [], often appear in string 
theories. Besides being of theoretical (academic, methodological) interest, the study on the 
thermodynamical aspects of them may be applied to the scenario of 'string cosmology'. 



II. 'BPS BLACK HOLES' IN SLOW MOTION 

We consider 'BPS black holes' in (d + 1) dimensional Einstein-Maxwell-dilaton theory, 
which is governed by the action 



S= d a+L x 



d+l„ V 9 



4 



4a 



R - ^—V^W - e-^F^F^ . (1) 



Here R is the scalar curvature and F^ = d^A v — d v A^ = 0,1, ... ,d), while is a 

dilaton field and a is a coupling constant. G denotes the (d + 1) dimensional Newton's 
constant. We set c = 1. 

This theory admits static multi-centered solutions whose metric and field configurations 
are given by |J 

ds 2 = -V'^^dt 2 + V^^dx 2 , (2) 

4a j. 2a 

e -—4> = yjz^ ? ( 3 ) 

^-Jw^b-v)' (4) 

A i (x) = (i = l,...,d), (5) 

in terms of a harmonic function 

t r/ \ . , 2{d-2 + a 2 ) 4tt _ m a 

V{x) = 1+ (d-l)(d-2) A^ G ^ RP' (6) 

where r a = x — x a , x a is the position vector of 'point particle' a and the labels of 'particles' 
a,b, . . . run over 1, . . . , N. A^-i = 27r d ^ 2 /Y{d/2) is the volume of a unit d — 1 sphere. 



3 In it is shown that they saturate a gravitational analogue of the Bogomol'nyi bound on their 
mass and charges. 
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The static point sources corresponding to the 'particles' can be described by the action 

1 = "E f ds * [m a e^ + e a A^—\ , (7) 

with the relation 



e a 2(d-2 + a*) 

m a ~i d-i • (8) 

This relation is just the extremity condition for spherically symmetric black holes in the 
Einstein-Maxwell-dilaton theory [ 25j26| . Thus the solution can be viewed as the one for an 



l N BPS black hole' system, where the electric force is cancelled by the forces mediated by 
the graviton and the dilaton. 

While there is no static force between the 'black holes', a velocity dependent force arises 
if one gives them small velocities v a . Many authors calculated the effective theory to 0(v 2 ) 
for various types of 'black holes'. One finds the energy of the 'BPS black holes' in the 
Einstein-Maxwell-dilaton theory [[12] 



H = E E v ak v bt (Sl5i + 5 M 5v - 5i6i)d ai d bj L , (9) 

a b 

where 

1 /* j 2(d-l) 

L = ~32^gJ dxV ^ {xh (10) 
with V given by @ and the spatial indices i, j, . . . run over 1, . . . , d. 

III. THE CANONICAL PARTITION FUNCTION 

The canonical partition function for N identical particles at temperature T in d spatial 
dimensions is 



4 If a dilaton field is coupled, the extreme limit of the black hole solution has a singularity in 
general. Therefore, we use the 'quotation marks' around the BPS black hole in this paper. 
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Zn = mw I [dq][dp] exp { ~ m ' 

where [dp] = Yla=i d d p a and [dq] = Yla=i d d x a . (3 = 1/T, where h is Planck's constant while 
Boltzmann's constant is set to be unity. 

Because there is no static forces, the Hamiltonian H for the system of iV 'BPS black 
holes' can be written as 

# = ^££ v A g AB v B = \v T gv, (12) 

A B Z 

where the label {A} = {(ai)} denotes the combination of the label of the particle and the 
spatial index. 

The canonical momentum can be found as 

arr 

Pa = %j [ = 'EGabv b = (Qv) a , (13) 

B 

and then the Hamiltonian is reexpressed as 



ff= ee^= yc-^. (M) 

A B 



Thus the partition function becomes 

2 /[dq]Vdet<7. (15) 



Nd 

y 1 /27T 



N\ \(3h 2 
If we define 



2(d-l) 4 47T — . 

yjz^i i x \ _i_ Q m y 

v ' A _ O A . . ^ 



d-2A d _ 1 V \r c \ d - 2 



(16) 



with 



v({C) = 1+ (d-i)(d-2) ^r Gro ? rp' (17) 

the Hamiltonian for N 'BPS black holes' with a common mass m can be rewritten as 

H = £™ £ + £ £ ^"(M ' + - SiSi)d m d bj L' . (18) 
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Then we may use 

G(ak)(be) = m 



SobSkt + -(8uS ij + $5} - 5{5\)d ai d hj L' 
m 



= m 



9(ak){bi) 



(19) 



and obtain the expression for the partition function 

Nd Nd 

1 (2Tim\ 2 /- ri . r 1 (2ixm\ 2 ^ 



(20) 



which is proportial to the volume of the moduli space Q. This feature is the same as that of 
point-particle systems in two (spatial) dimensions |27],[^|, in which there is no static force 
either. 

In the present system, the free energy F is expressed as 



F = -T\nZ = NT In iV — NT — —NT In — NT In V — T In -yj^ 

2 tj ft y/ 

Therefore the internal energy E is obtained as 

E 1 d((3F) d 



N N d/3 2 
which is the same as that of an ideal gas. 



T. 



(21) 



(22) 



IV. THE PERTURB ATIVE APPROACH 

To obtain the equation of state, we have to evaluate the volume of the moduli space Q. 
First, we try to obtain the volume of the moduli space perturbatively. One can rewrite 
this as 



where 



Q 



[dqj^detflf = J [dq] exp ^-Trlnflf 
[dq] exp -Tr(it 
H — Trn ■ 



[dq] 



-u 2 + 



-Tiu 2 + -(T™) 2 + 

4 8 V ; 



(23) 



uiakMbi) = —(5k£$ tj + 5 % k 5\ - 5 3 k 5\)d ai d hj L' . 
m 



(24) 



The term including the trace takes the form: 



2d 

Tru = — V d ai d ai L' 
m „ 



(25) 



and 



1W = — 



d W + 2(8 U 8" - <W ) E d ai d bj L'd ai ,d bjl L> . 



(26) 



a b 

To obtain the leading and next-to-leading contributions for a small m, we expand U as 

2 



L' = —- 



32tt 



1/ 4 4?r V ^ ~ ° 2 ^ 2 1 
2 [—2 AT; W^Tf™ SSuTui 



1 



2 j^-2 



+ 



1/4 4vr \ d - a 2 1 - a 



6 V^- 2A d _ x y 2(d- 1) d 
+0(m 4 ) . 



1 1 



c d e 



- 2 I f* ^ I ^ — ^ I I ^ — ^ 



(27) 



Using this, the traces can be written as f] 

,/ 4 47r \ d - a 2 
Tiu = d 



d- 2A d _ 1 I 2(d- 1) ' — v — ' |cc tt — a? 6 j 

a J= b 



d-2 



d ( 4 47T \ 2 d - a 2 1 - a 2 



+ - 



G 2 m 2 



2 Vd- 2 A d _iy 2(d- 1) d - 1 
1 1 



x 



1 



EEE _ ~Ad-2 \~. _™|d-2~'~EE| __ 2!</-2) 
/ / \' Xj a • Aj b\ \ Uj a "^c| / |*^a ^b\ 

b J= a,c f= a aj= b 



+0{{Gmf) 



(28) 



Tr u 2 = d 



4 47r\ 2 /d-a 2 \ 2 



G 2 m 2 



d-2A d _ 1 J \2(d-l) 
1 



x 



1 x x ^ 1 

+ EE 



EEE I _ ,/-2 I™ _ ™ I./-2 ,„ _ ™ 2U/-2) 

o ?= a,c a af=b 



+0{{Gm) 



(29) 



Then we can obtain the folowing expression for Q: 



5 See Appendix for detailed calculations. 
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Q_ 



d ( A An \ d - a 2 ^ , , 1 



d d X\d d x 2 

\X\ — X 2 \ d ~ 2 



d 



A 47T 



d - a 2 d - 2 + a 2 



A\d-2A d _ 1 J 2(d-l) 2(d-l) 



G 2 m 2 



x 



N{N-l){N-2) — 



d d x 1 d d x 2 d d x 3 



\xi — X 2 \ d ~ 2 \Xi — X 3 \ d ~ 2 



d 2 



X 



+2N(N - 1) 

( A 47T " 



1 

V 2 

2 



d — a 



d d X\d d x 2 
Xl -x 2 \ 2 ^~ 2 ) 

2 \2 

> G 2 m 2 



iV(iV-l)(iV-2)(iV-3) fi 



+4iV(iV-l)(iV-2) — 



d d x l d d x 2 \ 
X\ - x 2 \ d ~ 2 ) 
d d x 1 d d x 2 d d x 3 



\x\ — x 2 \ d ~ 2 \xi — x 3 \ d 2 



+2N(N-1)- 



d d x 1 d d x 2 
\ Xl - x 2 \ 2 ( d ~ 2 ) 



+0((Gm) 



(30) 



where V — J d d x. 

This expression includes divergent integrals and thus we must use the small cut-off length. 
These terms are, however, eliminated by 1/N factors in the N —>■ oo limit [[L]]. 
In the thermodynamical limit, or in the limit of large N, we obtain 
d( 4 4vr GmN\ d - a 2 i d ~ 2 ff d d x l d d x 2 



lim —In -^7 

N^oo N V N 



2 \d-2A d _ 1 i d - 2 ) 2[d-\) V 2 J J \xt 



x 2 



d-2 



d 2 



An GmN 



2 \d-2Ad-! £ d ~ 



+■ 



d 2 



An GmN 



V 



d — a 2 
.2Cd=l), 
d — a 2 



ffd-2 



V 2 



d d X\d d x 2 
\x\ — x 2 \ d ~ 2 



2 \d-2A d ^ t d ~ 2 ) \2(d-l) 
£ 2(d-2) rrr d d x 1 d d x 2 d d x 3 



d-2 + a 2 ' 
2d(d-a 2 ) 



x- 



V 3 



\ x l — x 2\ d 2 \ x l — x 3\ d ~ 2 



+- ^ i d - 2 ) ) 

where we have introduced a length scale I and the value N/£ d ~ 2 is fixed. 
For a spherical box of radius £, on can find 

ff d-2 ,, d d x 1 d d x 2 2d 



(31) 



V 2 



£Ci - X 2 



I d-2 



d + 2 



(32) 
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where V = A d _ x i A jd. 

The pressure P is derived from the free energy F as 

dF 



P 



dV 



thus we find 



(34) 



(35) 



Substituting (^TJ), (|3"2"D, (|3"3| ) into (EST), we obtain the equation of state in the limit of large 
N: 

PV 



NT 



with 



/( 2 /) = l-6 1 y + 6 2 y 2 + 0( 2 / 3 ) 



where 



61 



d-2 



b 2 = d{d - 2) 



~2{d-a 2 )' 


2 


ff 2d ) 


d- 1 




\d + 2J 



and 



2{d-a 2 )A 2d 
d - 1 J d + 2 ' 

f _ d-2 + a 2 ' 
[ ~ 2d{d - a 2 ) J (d + 2)(d + 4) 



An GmN 



d(5d + 2) 



(36) 

(37) 

(38) 
(39) 

(40) 



d-2A d _ 1 £ d ~ 2 ' 

is a dimensionless parameter, y can be regarded as a normalized coupling constant of the 
interaction. 

Now we can rewrite this equation of state in the Van der Waals form 

PV = NT with V = V [fiyT 1 , (41) 

which exhibits an 'anti-excluded' volume effect for a small y if a 2 < d. There is an effective 
attractive force between 'BPS black holes' for a 2 < d. The effective force is repulsive for 
a 2 > d, while for a 2 = d there is no interaction among the 'BPS black holes'. 
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V. THE STRONG COUPLING LIMIT AND PADE APPROXIMATION FOR THE 

EQUATION OF STATE 



In this section, we consider the strong coupling limit, y ^> 1. By counting the maximal 
number of G appearing in Q, we find 

Q Nd d-a 2 

— ~ y 2 x cons t, ; (42) 

since Q is an (Nd) x (Nd) matrix. [] 
This leads to, in the large y limit, 

PV d-2 d - a 2 _ da 2 - (d - 2) 2 

iVT^ 2~ d-2 + a 2 ~ 2(d-2 + a 2 ) ' ( ^ 

If a 2 < ^~p-, the pressure P becomes negative! This means that the gas will be sponta- 
neously condensed for large y. 

This feature can be explained by considering the moduli space metric for two body system 



II]. In the center-of-mass system, the moduli space metric in the small distance limit is 



ds^ f^^(dr 2 + rW), (44) 



where r is the distance between two particles and Q represents the angular variables. One 
can use a new radial coordinate near the origin: 

r (d-2) (d-a 2 ) 

R= drr i(d-2+a 2 ) . (45) 



The range of R is [0,oo) for a 2 > U=£L, while ( — oo, oo) for a 2 < d ■ Therefore in 
the case of a 2 < ^~ 2 ^ , two particles merges, or R — > — oo, for a sufficiently small impact 
parameter. This appears to be the cause of the instability of the gaseous state. 
An approximation which interpolates the small and large y regions shows 

O AT Nd d-a 2 

^ w (q P x(v)) = (1 + a iy )~^ , (46) 



6 Although the rank of Q may be less than Nd, this is a reasonable aproximation for N — > oo. 
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with 

2(d-2 + a 2 ) 2d 
^ = d-l d + 2- (47) 

Although this approximation does not utilize the information of the second order in 
the perturbative expansion, this form is very useful and provides the lowest-order Pade 
approximation for f(y) (see below) besides it garantees the positivity of the volume of the 
moduli space in the simplest way. 

In this approximation, the equation of state reads 

PV n _ (d-2)(d-a 2 ) g iy _ _ (d-2)(d-a 2 ) 2d y 

NT I[y) 2(d-2 + a 2 ) l + a lV d-l d + 21 + a lV ' 1 ' 

At the critical y = y c , y c satisfies f(y c ) = 0, the gas cannot remain in the normal gaseous 
phase. In the present approximation scheme, the value of y c is 

_ d + 2 d-l 
Vc ~ 2d (d-2) 2 -da 2 ' ( } 

for a 2 < a 2 c = (d - 2) 2 /d. 

The next-to-leading Pade approximation for the equation of state is also possible. The 
corresponding form for the moduli space volume is written by 



AT o 1 v u u,— u, 

(q p2 (y)f = (l + 2a iy + a 2 y 2 ) * , (50) 



with 



2 2 d-2 + a 2 , 

a * = 2al -—2^r^ h - (51) 



VI. THE MEAN-FIELD APPROXIMATION 



In |P-2|| , the effective field theory of 'BPS black holes' was constructed. Here we use the 
effective field theory to obtain the equation of state for a gas of 'BPS black holes' by the 
mean-field approximation. 
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In the mean-field approximation the number density with a chemical potential /i is given 



by 



n = — I d"p exp 



-P 



P 



,2mV d ~ 2 + a ' 2 

where the gauge interaction term is discarded. [] 



d{d-ar) 



(52) 



The mean-field potential V satisfies the following equation [[12 

d(d-a 2 ) 



d l V + 8tt 



d - 2 + a 2 
d- 1 



Gm n V 2 ( d - 2+aZ '> = 



(53) 



where we set n = uqV 2( - d - 2+a ' 2 ^ with a constant uq. For the spherical symmetric case, the 
equation for V reads 



d 2 V d-ldV 2d{d-2)(d-2 + a 2 ) Vn 



dr 2 r dr 



- + 



d-1 



V L_^.y2(d-2+a 2 ) — 

y N 



(54) 



with the boundary conditions 



dV 
df 



r=0 



dr 



2(d - 2)(d - 2 + a 2 
d-1 



(55) 



and 



V\ f=x = 1 + 



2(d- 2 + a 2 
d - 1 



-y 



(56) 



Here r = r/£ and we used / d d x n = N in the spherical box of which radius is i. 

To obtain the relation between y and N/(rioV), one can use the scaling covariance. It 
can be found that one has to solve 



d-1 , 2d(d - 2)(d - 2 + a 



d-1 



2\ d(d-a z ) 
t r2(d-2 + a 2 ) 

v 



(A) = 



(57) 



with the boundary conditions 



Vtf(O) = 0, V (0) = 1 , 



(58) 



7 Non-perturbative effect of the interaction was studied in [|T^]. 
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instead of solving (|54]) with ([55]) and (|56|). Here the prime (') denotes the derivative with 
respect to A. 

Next, one must solve the equation 

V>(\ s )_ (d-2) 2 -^±^y 



Xs v (\ s 

to get X s (y). 

Finally, N/ (noV) is given by 



1 , 2(rf-2+a 2 ) 

L ^~ d-i y 



(59) 



iV 



Qmf{y) 



d- 1 



n V 2(d - 2)(d - 2 + a 2 ) X 

In other words, solving we have 

d- 1 



V ° ,(As) f %(A.) + 



d-2 



A, 



(60) 



y 



(61) 



and this relation with (|60|) gives an implicit functional expression in terms of the parameter 
X s . 

Note that for a small y, we can find 



. , 2d 2 (d — a 2 ) o N 

w(8/) = n- ( ,_ l 1)( , + ' 2) y + o(y 2 ) 

which coincides with q p i(y) up to the first order in y. 



(62) 



d(d-q z ) 



On the other hand, for y grows as Vo(A s ) + x^^o'(^s) ~^ 0? we find Qmf(y) ^ y 2 ( d - 2+a ^ in 
the limit of large y. This behavior agrees with the previous analysis on the strong coupling 
region. 

In the mean-field approximation, the canonical partition function can be expressed as 



z _ 1 

Zn ~ m 



I ld d x 



\ - n N 



\(3h' 



- V d-2+a' 1 



1 / 2-Km 



[Qmf(y)] N ■ 



N\ \ I3h 2 



(63) 



Then, the equation of state reads PV/(NT) = f(y) with 

f(y) = 1 - d —^y^\nq mf (y). 



(64) 
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We evaluate f(y) by solving Eq. fl57|) by a numerical method, since the solution cannot 
be given in analytically in general. [| 

In Fig. [I], f(y) obtained by approximations considered above is shown. All the mean- field 
and Pade approximations indicate the consistent behavior of the function. The difference 
appears to be small when d — a 2 is small. 

There are regions where the pressure decreases with decreasing volume, where the system 
cannot be stable and will compress itself to a smaller volume. The isothermal compressibility 
K, 

K =-vw = M m+i ^ vm T' (65) 

becomes negative in such a region. 

The specific heat at constant pressure cp given by 

cp = c v - ^i|P- = \ + (f(y)f (f(y) + ^V(V)) 1 , (66) 

also becomes negative in such a region. 

These two quantities diverges when y — > y , and become negative for y > y . 

The values of y c (f(y c ) = 0) and y (K~ l = Cp 1 = 0) for several cases are listed in Table 1. 
(mf) means that the value is obtained by the mean-field aproximation, while (pi) and (p2) 
mean that the values are estimated by the leading and next-to-leading Pade approximations. 



? See, however, for analytic approximations, |29,30|. See also [31] 
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d = 3 


d = 4 


d = 5 


d = 5 




a 2 = 


a 2 = 


a 2 = 


a 2 = l 


Vc (mf) 


1.33354 


0.374997 


0.17789 


0.399981 


Vc (p2) 


1.47578 


0.46414 


0.245029 


0.549095 


Vc (P 1 ) 


1.66667 


0.562500 


0.311111 


0.700000 


Vo (mf) 


0.944025 


0.248672 


0.120393 


0.235276 


yo (p2) 


1.04266 


0.288196 


0.144323 


0.292631 


l/o (pi) 


1.17851 


0.347644 


0.180799 


0.374393 



Table 1 



We have obtained the similar critical values for y by the different approxmations. Thus 
the precise values are expected to be in the vicinity of the approximated values. This will 
be confirmed by the numerical simulation in the future. 

VII. PARTICLE DISTRIBUTIONS 

We consider a spherical distribution of the 'BPS black holes'. We define the mass distri- 
bution in the mean-field approximation by 

M(f) = A d _! ^ \V(x)}^-Za2) X d-^dx . (67) 
Jo N 

M(f) gives the fraction of the mass inside the sphere with the radius f. 
This can be approximated as 

M(f) w r D , (68) 

noting that M(0) = and M(l) = I. 

D is considered as the most naive definition of the fractal dimension. We evaluate the 
value of D at the least mean square of the deviations. The results for d — 3, 4, 5 are exhibited 
in Table 2, 3 and 4. 
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d = 3 


y 


0.2 


0.4 


0.6 


0.8 


1.0 


1.2 


1.4 


1.6 


1.8 


2.0 


D (a 2 = 0) 


2.69 


2.46 


2.29 


2.15 


2.04 


1.95 


1.87 


1.81 


1.75 


1.70 


D (a 2 = 1/3) 


2.74 


2.58 


2.46 


2.38 


2.31 


2.26 


2.22 


2.19 


2.16 


2.14 


D (a 2 = 1) 


2.83 


2.74 


2.69 


2.65 


2.62 


2.60 


2.59 


2.58 


2.57 


2.56 



Table. 2 



d = 4 


y 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


D (a 2 = 0) 


3.52 


3.11 


2.76 


2.42 


2.09 












D (a 2 = 1) 


3.68 


3.45 


3.29 


3.17 


3.07 


3.00 


2.93 


2.88 


2.84 


2.80 



Table. 3 



d = 5 


y 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


D (a 2 = 0) 


4.00 


2.90 


















D {a 2 = 1) 


4.29 


3.77 


3.35 


3.02 


2.78 












D (a 2 = 9/5) 


4.48 


4.14 


3.91 


3.75 


3.62 


3.53 


3.45 


3.39 


3.33 


3.29 



Table. 4 



We find that D slowly decreases as y grows. These behaviors must be checked by the 
numerical simulations in the future. 

VIII. SUMMARY 

In summary, we have studied many-body system of 'BPS black holes' in d dimensions. 
The canonical partition function of the system is proportional to the volume of the moduli 
space of iV 'BPS black holes'. Therefore the difference from an ideal gas arises through the 
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change of the effective volume of the system: the internal energy and heat capacity at a 
fixed volume are the same as those of an ideal gas. 

We have estimated the equation of state of the gas by perturbations in the case of a 
weakly self-interacting gas. The appropriate variable for the expansion has been found to 
be y oc Gm/£ d ~ 2 , where £ is the length scale of the box containing the gas. We have also 
estimated the equation of state of the gas by Pade approximation associated with counting 
the highest order of the couplings, and by the mean-field method. 

In these approximation, we have found that the pressure simply decreases as the value 
of y increases if the temperature and the density are fixed. This feature is common in the 
cases for the different spatial dimensions and the dilaton couplings. 

For each dimensionality d, the critical value for the dilaton coupling a has been found. 
If the absolute value of the coupling is smaller than the critical value a c = ^2, there exists 
the critical value y c ; for y > y c , the pressure becomes negative. 

A usual thermodynamical limit, iV — > oo, V — > oo and N/V fixed, can be taken only if 
a 2 > a 2 . This limit can be obtained by taking y — > oo. The equation of state is then 

PV da 2 -(d- 2f 



NT 2{d-2 + a 2 ) ' 
The isothermal compressibility K is 



(69) 



K = ^d-2 + a 2 ) 

NT da 2 - (d - 2) 2 ' 1 ' 

and the specific heat at constant pressure cp is 

d da 2 -(d-2f 
Cf= 2 + 2( rf -2 + aV (71) 

Of course the analytical attempt to describe thermodynamic property of the 'BPS black 
hole' gas may have been inconclusive and numerical simulations of the self-interacting system 
may be needed. The present work provides a guideline for such numerical calculations to 
survey the critical points of the 'BPS black hole' gas. 
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For the near-extremal cases, we may treat static potentials as a perturbation from the 
critical coupling. These cases will be interesting because they may exhibit complicated 
temperature-dependent critical behaviors. 

ACKNOWLEDGEMENTS 

KS would like to thank Kenji Sakamoto for reading this manuscript. He also thanks 
Yoshinori Cho for useful comments. 



18 



APPENDIX: SOME USEFUL FORMULAS 
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FIG. 1. f(y) is plotted against y for d = 3, d = 4, d = 5 and d = 9. The solid line is 
obtained by the mean-field approximation, the broken line by the lowest-order Pade approximation, 
the dot-broken line by the next-to-leading Pade approximation, and the dotted line is the series 
expansion up to the second order. 
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